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Abstract 

In this paper, Picard method is proposed to solve the Cauchy reaction-diffusion equation with fuzzy initial condition 
under generalized /-/-differentiability. The existence and uniqueness of the solution and convergence of the proposed 
method are proved in details. Some examples are investigated to verify convergence results and to illustrate the 
efficiently of the method. Also, we obtain the switching points in examples. 
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Introduction 

As we know the fuzzy differential equations FDE are 
one of the important part of the fuzzy analysis theory 
that play major role in numerical analysis. For example, 
population models (Guo et al. 2003), the golden mean 
(Datta 2003), quantum optics and gravity (El Naschie 
2005), control chaotic systems (Feng and Chen 2005; Jiang 
2005), medicine (Abbod et al. 2001; Barro and Marin, 
2002). Recently, some mathematicians have studied FDE 
(Abbasbandy and Allahviranloo 2000; Abbasbandy et al. 
2004; Abbasbandy et al. 2005; Allahviranloo et al. 2007; 
Bede 2008; Bede and Gal 2005; Bede et al. 2007; Buckley 
and Feuring 2000; Buckley and Jowers 2006; Buckley et 
al. 2002; Chalco-Cano and Romn-Flores 2006; Chalco- 
Cano and Romn-Flores et al. 2007; Chapko and Johansson 
2012; Chen and Ho 1999; Cho and Lan 2007; Congxin 
and Shiji 1993; Diamond 1999; Diamond 2002; Ding et 
al. 1997; Dubois 1982; Dou and Hon 2012; Fard 2009a,b; 
Fard and Bidgoli 2010; Fard and Kamyad 2011; Fard et al. 
2009; Fei 2007; Jang et al. 2000; Jowers et al. 2007; Kaleva 
1987,1990,2006; Lopez 2008; Ma et al. 1999; Mizukoshi 
et al. 2007; Oberguggenberger and Pittschmann 1999; 
Papaschinopoulos 2007; Puri and Ralescu 1983; Seikkala 
1987; Solaymani Fard and Ghal-Eh 2011; Song et al. 2000). 
The fuzzy partial differential equations FPDE are very 
important in mathematical models of physical, chemical, 
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biological, economics and other fields. Some mathemati- 
cians have studied solution of FPDE by numerical meth- 
ods (Afshar Kermani and Saburi 2007; Allahviranloo 2002; 
Barkhordari Ahmadi and Kiani 2011; Buckley and Feuring 
1999; Chen et al. 2009; Farajzadeh et al. 2010; Moghadam 
and Jalal 2011; Rouhparvar et al. 2010; Verma et al. 2009). 
In this work, we present the Picard method to solve the 
Cauchy reaction-diffusion equation as follows: 

Ut (x, t) = u xx (x, t) + u(x, t), 0 < t < T, 

a < x < b, a,b,T e R. 

(1) 



With fuzzy initial condition: 
u(x, 0) = f(x). 



(2) 



Sprin 



The structure of this paper is organized as follows: In 
section "Basic concepts", some basic notations and defini- 
tions in fuzzy calculus are brought. In section "Descrip- 
tion of the method" Eqs.(l,2) are solved by Picard method. 
The existence and uniqueness of the solution and con- 
vergence of the proposed method are proved in section 
"Existence and convergence analysis" respectively. Finally, 
in section "Numerical examples", the accuracy of method 
by solving some numerical examples are illustrated and a 
brief conclusion is given in section "Conclusion". 

Basic concepts 

Here basic definitions of a fuzzy number are given as 
follows, (Allahviramloo 2005; Dubois and Prade 2005; 
Kauffman and Gupta 1991; Nguyen 1978; Zadeh 1965) 



ger 
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Figure 1 The results of Example 5.1 for (m(x, 0.6, 0.1), m(x, 0.6, 0.1)). 



Definition 2.1. An arbitrary fuzzy number 5 in the para- 
metric form is represented by an ordered pair of functions 
(u, u) which satisfy the following requirements: 

(i) u : r — > u(r) e R is a bounded left-continuous 
non-decreasing function over [0,1], 

(ii) u:r^- u(r) e R is a bounded left-continuous 
non-increasing function over [ 0, 1], 

(iii) u(r) < u(r), 0 < r < 1. 

Definition 2.2. For arbitrary fuzzy numbers u, v e E , 
we use the distance (Hausdorff metric) (Goetschel and 
Voxman 1986) D(u(r), v(r)) = max{sup re [ 01 ] \u(r) — 
v(r)|, sup \u(r) — v(r)\}, and it is shown (Puri and Ralescu 
1986) that (E 1 , D) is a complete metric space and the 
following properties are well known: 



D(u + w,v+ w) = D(u,v),Vu,v g E 1 , 

D(ku, kv) =| k | D(u,v), V k e R, u,v e £\ 

D(u + v, w + e) < D(w, w) + D(v,e), V w,v, w,e e f 1 . 



Definition 2.3. Consider 3c,j/e £. If there exists z e E 
such that jc = j + z then z is called the /f- difference of X 
and y, and is denoted by x 0 y (Bede and Gal 2005). 

Proposition 1. If / : {a, b) — > E is a continuous fuzzy- 
valued function then g(x) = f*f(t) dt is differentiable, 
with derivative g (x) = fix) (Bede and Gal 2005). 



Definition 2.4. (see (Bede and Gal 2005)) Let/ : {a, b) -> 
E and xq € (a, We say that / is generalized differen- 
tiable at xq (Bede-Gal differentiability), if there exists an 
element / (#rj) e E, such that: 

i) for all h > 0 sufficiently small, 3 f(xo + h) Q f(xo), 
3 f(xo) 0 f(xo — h) and the following limits hold: 

,. f(x 0 + h)ef(x 0 ) f(x 0 )Qf(x 0 -h) 

lim = lim 

h^0 h h^O h 

= f (xo) 

or 

ii) for all h > 0 sufficiently small, 3 f(xo) 0 f(xo + h), 
3 f(xo — h) 0 /0*o) and the following limits hold: 

/ (xo) 0 f(x 0 + /z) _ /(jcq -h)Q f(x 0 ) 
—h h^o —h 

= f (xo) 



lim 

h^0 



or 



iii) for all h > 0 sufficiently small, 3 /(#o + h) 0 f(xo), 
3 /C*o — h) Q f(xo) and the following limits hold: 



lim 



/(*o + fe) 0 /(jcq) 



lim 

h^0 



f(x 0 - h) ef(x 0 ) 



-h 



=f (x 0 ) 



or 



iv) for all h > 0 sufficiently small, 3 f(xo) 0 f(xo + /z), 
3 /(^o) © f(Xo ~ h) an d the following limits hold: 



lim 

h^0 



f(xo) 0 f(xo + h) 



lim 

h->0 



f(xo) ef(x 0 - h) 



--f (xo) 



Definition 2.5. Let / : {a, b) E. We say / is (i)- 
differentiable on (a, b) if / is differentiable in the sense 
(i) of Definition (2.7) and similarly for (ii), (iii) and (iv) 
differentiability. 

Definition 2.6. (see (Chalco-Cano and Romn-Flores 
2006)) Let / : (a,b) ->■ E and #o G («, b). A point 
G (a, b) is^said to be a switching point for the differ- 
entiability of /, if in any neighborhood V of xo there exist 
points x\ < Xq < X2 such that 

(type 1)/ is differentiable at x\ in the sense (z) of 
Definition (2.6) while it is not differentiable in the 
sense (ii) of Definition (2.6), and / is differentiable at 
X2 in the sense (ii) of Definition (2.6) while it is not 
differentiable in the sense (i) or Definition (2.6), or 

(type 2) / is differentiable at x\ in the sense (ii) of 
Definition (2.6) while it is not differentiable in the 
sense (i) of Definition (2.6), and / is differentiable at 
X2 in the sense (i) of Definition (2.6) while it is not 
differentiable in the sense (ii) or Definition (2.6). 
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Proposition 2. (see (Chalco-Cano and Romn-Flores 
2006)) Let/ : (a, b) -+ E and x 0 e (a, b). 

(a) If Xo e (a, b) is a switching point for the 
differentiability of / of type 1, then / is differentiable 
at xq in the form (iv). 

(b) If Xq € (a, b) is a switching point fo.r the 
differentiability of / of type 2, then / is differentiable 
at xo in the form {Hi). 

Definition 2.7. A triangular fuzzy number is defined as a 
fuzzy set in E , that is specified by an ordered triple u = 
(a,b,c) e R 5 with a < b < c such that u(r) =[u(r),u(r)] 
are the endpoints of r-level sets for all r €[0,1], where 
u(r) = a + (b — a)r and u(r) = c — (c — b)r. Here, w(0) = 
a,u(0) = c,u(l) = 5(1) = b, which is denoted by u(l). 
The set of triangular fuzzy numbers will be denoted by E l . 

Definition 2.8. (see^ (Chalco-Cano and Romn-Flores 
2006)) The mapping / : T — > E" for some interval T 
is called a fuzzy process. Therefore, its r-level set can be 
written as follows: 

(f(t))(r)=[f(t,r),f(t,r)], t e T, re[0,l]. 

Definition 2.9. (see (Chalco-Cano and Romn-Flores 
2006)) Let / : T -> E" be Hukuhara differentiable 
and denote (f(t))(r) =[f(t, r),f(t, r)]. Then, the boundary 

function f(t, r) and f(t, r) are differentiable (or Seikkala 
differentiable) and 

(f'(t))(r)=[f(t,r),f(t,r)], t e T, re[0,l]. 
If / is (ii) -differentiable then 

(f'(t))(r)=[f (t,r),f(t,r)], tef, re[0,l]. 

Description of the method 

To obtain the approximation solution of Eqs.(l,2), based on 
Definition (2.6) we have two cases as follows: 

Case (1): u(x, i) is (z) -differentiable, in this case we 
have, 

u(x, t) = f{x) + / [u(x, t) + u xx (x, t)\ dt. (3) 
Jo 

Case (2): u(x, i) is («') -differentiable, in this case we 
have, 

u(x, t) =f(x) e (-1). / [ u(x, t) + u xx (x, t)] dt. 
Jo 

(4) 

Now, we can write successive iterations (by using Picard 
method) as follows: 



Case (1): 
uo(x, t) =f(x), 

u n+ i(x, t) =f(x) + Jq[u„(x, t) + u„ xx (x, t)] dt, n > 0. 

(5) 

Case (2): 
u 0 (x, t) =f(x), 

u„+i(x, t) =f(x) e (-1). / [u„(x, t) + Un xx (x> dt, 
Jo 

n > 0. 
(6) 

Remark 1. For u xx we have cases as follows: 



Case (1): u and u be (i)-differentiable and u and u be 
(ii) -differentiable 

u xx {x, t) =[ Ux X (x, t, r), u xx (x, t,r)]. 

Case (2): u is (z)-differentiable and u is 

(ii) -differentiable and u is (w)-differentiable and u is 

(i)-differentiable 

u" (x,t) =[ u" (x, t, r) , u (x,t,r)]. 

Remark 2. We discuss about switching points as follows: 

Case (1): u is (z')-differentiable 

If 3^ <0) ^^ >0,x€[a,xo] and 

> o, M^Tl < o, *e(«o,*]andre[ai] 
then we have u (x, t) =[ u (x, t, r), u (x, t, r)] and Xo is 
a switching point in the form (iv). 

<0, 3 -^H >o,*e[«,*i]and 



3a: 2 ' 3x 2 

d 2 u(x,t,r) r. d 2 u(x,t,r) 

> o, - 



-^i— - ». dx 2 < 0, x e (xi,b] andr e[0,l] 
then we have u (x, t) =[ u (x, t, r), u (x, t, r)] and X\ 
is a switching point in the form (iv). 
Case (2): u is (w)-differentiable 

, c du(x,t,r) „ du(x,t,r) n r l j 
If dx > °' Sx < °' X € \- a > x 01 and 
3^ < a Suix^r) >0) xe (xQib] and r e[0> x] 

then we have w' f) =[ «' (x, t, r), u (x, t, r)] and xo is 
a switching point in the form (Hi). 



If 



dx 2 



> 0, 



dx 2 



< o, 



< 0,x €[a,Xi] and 
1^2— - », . )x2 > 0, * e andr e[ 0,1] 
thenwehavew (x,t) =[u (x,t,r),u (x, t, r)] andxi 
is a switching point in the form (Hi). 
Case (3): u is (z')-differentiable 
If aufgr) < 0j 3»^r) > o, y x e [ a,b] then we have 



dx 



u (x, t) =[u (x, t, r), u (x, t, r)]. 



If 



3* 2 



< 0, 3 E g*2' r> > 0,x e[a, xq] and 



d 2 u(x,t,r) 
dx 2 



0, 3 - 3 ( ^ t,r) < 0, *e (*o,&] andr e[ 0,1] 
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then we have u (x, t) = [ u" (x, t, r), u (x, t, r)] and xq 
is a switching point in the form (iv). 
Case (4): u is («7)-differentiable 
If^f >0, d -^ <0,V*e[ a,b] then we have 
u (x, t) =[ u (x, t, r), u (x, t, r)]. 
If%^>0,%^<0,*e[«,*i]and 
%^<0,%^>0, and re[0,l] 

then we have 2 (#,£)=[« (*> t,r),u (x, t, r)] and #i 
is a switching point in the form (Hi). 

Existence and convergence analysis 

In this section we are going to prove the existence and 
uniqueness of the solution and convergence of the method 
by using thefollowing assumptions. 
Consider f(x) is bounded for all x € [ a, b] and 



Table 1 Numerical results for Example 5.1 



D(u xx (x, t), u*(x, t)) < LD(u(x, t), u* (x, t)), L > 0. 



Let, 



a = (T + TL). 

Lemma 1. \iu,v,w e E n and X e R, then, 

(i) D(uQ v, uQ w) = D(v, w), 

(ii) D(QXu, QXv) =| 1 | D(u,v). 

Proof (i). By the definition of D, we have, 
D(uQv,uQ w) 

= max{sup rG[01] | u(r) - v(r) - u{r) - w(r) |, 

su Pre[o,i] I »W - Kr) -u(r)- w(r) |} 
= max{sup,. G[(U] | (u(r) - v(r)) - (u(r) - w{r) \, 

su Pre[o,i] I - v(r)) - (u(r) - w(r)) |} 
= max{sup re[ai] | w{r) - v(r) |,sup rG[01] | w(r) - v{r) |} 
= max{sup rG[01] | v(r) - w(r) |,sup rG[01] | v{r) - w(r) |} 
= D(y,w). □ 

Proof (ii): 

D{ 0 ku, QXv) 

= max{sup rG[011 | lM(r) - A.v(r) |, 



sup re[0 ,i]|A.M(r) - Xv{r) |} 
= max{sup rG[0;1] | Xu(r) - Xv(r) \, 

su Pre[0,l] I Xu ( r ) ~ Mr) |} 
= D(Xu,Xv) =| X | D{u,v). 



X 


(u,r = 0.1, n = 25, t = 0.6) 


(u,r = 0.1, n = 25, f = 0.6) 


-0.2 


0.34245 


0.68532 


-0.1 


0.35571 


0.67259 


0.0 


0.36438 


0.66135 


0.1 


0.35724 


0.67056 


0.2 


0.34862 


0.67843 


0.3 


0.33597 


0.68453 


0.4 


0.32551 


0.69064 



Theorem 1. Let 0 < a < 1, then Eqs.(l,2), have 
an unique solution and the solution u n (x, i) obtained 
from the relation (8) using Picard method converges to 
the exact solution of the problems (1,2) when u is (ii)- 
differentiable. 

Proof. Let u and u* be two different solutions of Eqs.(l,2) 
then 

D(u(x, t), u*(x, £)) 

= D(f(x) Q (-1). f [u(x, t) + Uxx(x, t)} dt, 
Jo 

f(x) 9 (-1). / [u*(x, t) + u* xx (x, t)} dt) 
Jo 

= D(e(-l). [u(x,t) + u xx (x,t)] dt, 
Jo 

0(-l). / [u*(x,t)+u* xx (x,t)] dt) 
Jo 

< T(D(u(x, t), u* (x, t))) + TL(D(u(x, t), u* (x, t))) 
= aD(u(x, t), u*(x, t)). 

From which we get (1 — a)D(u(x, t), u* (x, t)) < 0. Since 
0 < a < 1, then D(u(x, t),u*(x, t)) = 0. Implies u(x, t) = 
u*(x, t). 

Also, we have 

D(u n+ i(x, t),u(x, t)) < aD(u n ,u). 

Since, 0 < a < 1, thenD(«„(^, t),u(x, t)) — >• 0 as n — >• oo. 
Therefore, u n (x, t) —*■ u(x, t). □ 

Table 2 Numerical results for Example 5.2 



□ 



X 


(o,r = 0.2,n = 17,f = 0.7) 


(i7, r = 0.2, n = 1 7, t = 0.7) 


0.1 


0.3726508 


0.7023467 


0.2 


0.3827526 


0.693941 1 


0.3 


0.3957162 


0.6848745 


0.4 


0.4136805 


0.6612719 


0.5 


0.4372286 


0.6559328 


0.6 


0.4438574 


0.6337306 
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Table 3 Numerical results for Example 5.2 



X 


(u,r = 0.2,n = 21, f = 0.7) 


<y,r = 0.2,n = 21, f = 0.7) 


0.1 


0.4646316 


0.7954225 


0.2 


0.4717726 


0.7844835 


0.3 


0.4823549 


0.7737607 


0.4 


0.5074658 


0.7528153 


0.5 


0.5263437 


0.7474326 


0.6 


0.5309183 


0.7215178 



Remark 3. The proof of other case is similar to the 
previous theorems. 

Numerical examples 

In this section, we solve the Cauchy reaction-diffusion 
equation by using the Picard method. The program has 
been provided with Mathematica 6 according to the fol- 
lowing algorithm where £ is a given positive value. 

Algorithm : 

Step 1. Set n <- 0. 

Step 2. Calculate the recursive relations (7) or (8). 

Step 3. lfD(u n+ i(x, t),u„(x,t)) < e then go to step 4, 

else n <— n + 1 and go to step 2. 

Step 4. Print u„(x, t) as the approximate of the exact 

solution. 

Example 5.1. Consider the Cauchy reaction-diffusion 
equation as follows: 

u t (x, t) = u xx (x, t) + u(x, t). (7) 

With initial condition: 

u(x, 0) =/(*) = ((1 - r)x 3 , (r - l)* 3 ). (8) 

e = 10 -4 . x = 0 is a switching point. 

Case (1): n = 22 and a = 0.8652. 
Case (2): a = 0.84569. 

Table 1 shows that, the approximation solution of 
the Cauchy reaction-diffusion equation is convergent 
with 25 iterations by using the Picard method when 
u is (ii)-differentiable. 

Example 5.2. Consider the Cauchy reaction-diffusion 
equation as follows: 

u t (x, t) = u xx (x,t) + u(x,t). (9) 

With initial condition: 

u(x,0) =f(x) = (x 2 + l,x 2 + 2,x 2 + 3). (10) 

€ = 10" 3 . 

Case (1): a = 0.7546. 

Table 2 shows that, the approximation solution of 
the Cauchy reaction-diffusion equation is convergent 



with 17 iterations by using the Picard method when 
u is (i)-differentiable. 
Case (2): a = 0.7762. 

Table 3 shows that, the approximation solution of 
the Cauchy reaction-diffusion equation is convergent 
with 21 iterations by using the Picard method when 
u is (ii)-differentiable. 

Conclusion 

The Picard method has been shown to solve effectively, 
easily and accurately a large class of nonlinear problems 
with the approximations which convergent are rapidly 
to exact solutions. In this work, the Picard method has 
been successfully employed to obtain the approximate 
solution of the Cauchy reaction-diffusion equation under 
generalized //-differentiability. 
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